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I We report the discovery of a 2; = 2.65 low-ionization iron broad absorption Hne 

quasar, FIRST J100424. 9+122922, which is gravitationally-lensed by a galaxy at 2; ~ 
ly-^ ■ 0.95. The object was discovered as part of a program to find very red quasars by match- 

ing the FIRST radio survey with the 2-MASS near-infrared survey. J100424.9+122922 
^ . is the second lensed system to be found in this program, suggesting that many gravi- 

^2 ■ tational lenses are probably missed from conventional optical quasar surveys. We have 

i made a simple lens model and a rough estimate of the reddening in the immediate en- 

I vironment of the quasar which suggests that the quasar is intrinsically very luminous 

I and is accreting at close to the Eddington limit of its ~ 10^ Mq black hole. The lens- 

' ing galaxy has a small amount of dust which is responsible for some excess reddening 

observed in the fainter image of the quasar, but is otherwise a fairly typical massive 
elliptical galaxy. We model the selection effects working against the detection of red 
quasars in both lensed and unlensed samples. We show that these selection effects are 



, very effective at removing even lightly-reddened high redshift quasars from magnitude- 

^ ' limited samples, whether they are lensed or not. This suggests that the red quasar 

population in general could be very large, and in particular the class of iron broad ab- 
^ ' sorption line quasars of which J100424. 9+122922 is a member may be much larger than 

• their rarity in magnitude- limited samples would suggest. 

Subject headings: quasars: individual (J100424. 9+122922) - quasars: absorption lines 
- gravitational lensing - radio continuum: galaxies 
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1. Introduction 

We are currently conducting a survey for very red objects selected using the Two-micron All 
Sky Survey (2MASS; Kleinmann et al. 1994) and the Faint Images of the Radio Sky at Twenty- 
centimeters survey (FIRST; Becker, White & Helfand 1995; White et al. 1997 ). We select stellar 
infrared objects coincident with FIRST radio sources which lack (or have only faint) counterparts 
on the Palomar Observatory Sky Survey plates. To date wc have examined about 100 candidate 
red quasars, most using Keck or Lick spectroscopy. Of these, ~ 13 have optical spectra charac- 
teristic of z > 0.5 quasars reddened by dust (final classification will need to await near-infrared 
spectroscopy as broad lines sometimes appear in the near-infrared where the extinction to the 
quasar broad line region is lower). In a previous paper (Gregg et al. 2002) we reported the discov- 
ery of two dust-reddened quasars discovered as part of this survey FIRST J013435. 6-093102 and 
FIRST J073820.1-F275045 - one of which (J0134-0931) is gravitationally lensed (see also Winn et 
al. [2002]). In this paper we describe the recent discovery of a second gravitationally- lensed system, 
FIRST J100424.9-M22922 (hereafter J1004-M229). This quasar is reddened by a combination of 
dust and extreme, low ionization broad absorption in the rest-frame ultraviolet, making it a mem- 
ber of the rare sub-class of low-ionization broad absorption line quasars, dubbed "FeLoBALs" by 
Becker et al. (1997), as much of the absorption is due strong absorption from metastable Fe II 
(Hazard et al. 1987; Becker et al. 1997). 

There arc two main hypotheses concerning the nature of broad absorption line quasars (BALs) . 
One is that most quasars possess BAL flows, but wc only see them if the quasars are viewed at 
a particular angle to the line of sight which grazes a torus of gas around the nucleus (hereafter 
the "orientation hypothesis", e.g. Weymann et al. 1991). The other, which applies particularly 
for the case of the low-ionisation BALs [LoBALs], is that BALs are quasars at an early stage in 
their evolution (hereafter the "youth hypothesis", e.g. Voit, Weymann & Korista 1993). Recently, 
support for the youth hypothesis has come both from an analysis of the radio properties of BALs in 
FIRST (Becker et al. 2000), which indicates a range of orientations can produce BAL quasars, and 
also an apparent association of BALs with infrared-luminous quasars having hosts with merging 
morphologies and host spectra with "post-starburst" features (Canalizo & Stockton 2001). A 
further, but more controversial, hint that BALs may represent a stage in the life cycle of a quasar 
is that BALs may have steep intrinsic X-ray spectral slopes (Mathur et al. 2001; but see also 
Gallagher et al. 2001 who argue that the X-ray spectral slopes of BALs are not significantly steeper 
than those of normal quasars). Steep X-ray slopes may be connected with high accretion rates 
relative to the Eddington limit (Laor et al. 1997), so if BALs do indeed have steeper X-ray slopes 
this would suggest a link between high accretion rates and the BAL phenomenon. If it is assumed 
that high accretion rates are typical in the early life of a quasar, perhaps because the supply of fuel 
capable of falling into the black hole is greatest just after the triggering event, then a consistent 
picture emerges in which BAL flows are a symptom of high accretion rates produced as radiation 
pressure from a young quasar close to the Eddington limit pushes away excess material. 

The first FeLoBAL to be identified was LBQS 0059-2735 (Hazard et al. 1987). Since then five 
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more have been found in the FIRST Bright Quasar Survey (FBQS) (Becker et al. 1997; 2000). These 

objects arc overlooked in quasar surveys based on UV/optical color or objective prism selection, 
not only because of their extreme absorption features and lack of prominent emission lines, but also 
because, like other LoBALs, they tend to be more heavily reddened by dust than the normal quasar 
population (Spraybcrry & Foltz 1992). Like normal BALs, few are very radio-loud, so they do not 
appear in surveys with high radio flux density limits either. Quasar surveys based on the FIRST 
radio survey, are, however, better suited to finding such objects in significant numbers, because of 
the faint radio flux limit and the lack of the necessity for candidate selection based on exceptionally 
blue colors. 

For our other FIRST/2MASS gravitational lens, J0134-0931, we were unable to determine 
whether the reddening arose predominately in the immediate environment of the quasar, or in 
the lensing galaxy. Determining the source of the reddening (and extinction) is critical, however: 
if heavy reddening by lenses is commonplace, then many gravitationally lensed systems may be 
missing from optically-based lens surveys. If, on the other hand, the reddening is in the host, then 
many red quasars may be missing from current quasar surveys, with important implications for 
attempts to match the mass accreted onto black holes during the "quasar epoch" to the masses 
of black holes in galaxies today (e.g. Merritt &: Ferrarese 2001). An association between BAL 
quasars and gravitational lensing was noted by Chartas (2000), suggesting that reddening in the 
hosts usually dominates, but examples of lensed quasars reddened by their lenses are also known 
[B0218+357 (e.g. Menten & Reid 1996), and PKS 1830-211 (e.g. Wiklind & Combes 1996)]. 

Our infrared imaging and optical and infrared spectroscopy demonstrate conclusively that 
J1004-I-1229 is a lensed FeLoBAL and allow us to derive the basic quantities of the lensing sys- 
tem and also investigate the nature of the lensing galaxy and quasar, as detailed in Sections 2-4. 
J1004-I-1229 is the second gravitationally lensed, z > 0.5 red quasar found in our survey. This 
lensing rate (2/13) is much larger than higher that for the overall quasar population, though the 
statistics are poor so far. In Section 5 we model the effects of lensing and reddening on the optical 
and near-infrared selection of quasars, attempting to account for this apparently high lensing rate 
of reddened quasars. 

We assume a cosmology with JIm = 0.3, J^a = 0.7, Hq = 60kms~^Mpc^^, unless otherwise 
stated. 

2. Observations 

Matching of FIRST and 2MASS showed that the K' = 14.5 2MASS source corresponding to 
J1004-I-1229 is 0.3 arcsec from the 12mJy FIRST radio source (compared to estimated rms position 
errors of 0.2 arscec in 2-MASS and 0.5 arcsec in FIRST), but that this object was not present on 
second-generation digitised sky survey plates to a limit of i? ~ 20.8. In this Section, we describe 
observations of J1004+1229 in the radio, near-infrared and optical bands. Broad band magnitudes 
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and radio flux densities for the system are listed in Table 1. 

2.1. Spectroscopy 

J1004+1229 was initially identified as an FeLoBAL in a Keck ESI spectrum taken on 2000 
December 29, though the gravitationally lensed counterimage was not noticed at the time. 

Infrared spectroscopy with the NASA Infrared Telescope Facility (IRTF) using NSFCAM was 
taken on 2001 April 16 (UT) using the ii'-band grism. The spectrum was nodded up and down 
the slit in the standard two-point "ABBA" pattern (Eales & Rawlings 1993) to facilitate sky 
subtraction. A iiT-band acquisition image showed two components separated by 1.5-arcsec which 
were fortuitously almost aligned with the north-south slit. The spectra show that both components 
have an emission line at an identical redshift (Figure 1). We denote the bright image as component 
A, and the faint image as component B throughout this paper. 

A further ESI spectrum was taken on 2001 May 26 with the 1-arcsec slit rotated to a sky PA 
of 9 deg. to include both images of the quasar (Figure 2). Two 900s integrations were taken, and 
both components of the lensed system were easily resolved along the slit in the QfIS arcsec seeing. 
The data were calibrated using an observation of Fiege 34. 

Figure 3 compares the ESI spectrum of J1004+1229 with two other FcLoBALs, LBQS 0059- 
2735 (Hazard et al. 1987) and FIRST 1556-^3517 (Becker et al. 1997). The similarity of J1004-M229 
and FIRST 1556-f3517 is evident. Further discussion of the spectrum may be found in Section 3.1. 
The spectral resolution of the observations is ~ 2A, though the spectra have been smoothed to 
8 A resolution in Figures 2 and 3. 

2.2. Optical and near-infrared imaging 

Infrared images in J and K' were taken with NSFCAM at the IRTF on 2001 April 16. The data 
were taken in a 9-point mosaic pattern with 60s exposures at each point, and were reduced using 
the DIMSUM package after flat-fielding with dome flats. The weather was non-photometric and the 
seeing 0'.'7. V— and /—band images were obtained with the Andalucia Faint Object Spectrograph 
Camera (ALFOSC) on the Nordic Optical Telescope (NOT) on 2001 April 22 for 3 x 300s and 
3 X 600s respectively, and flat-fielded with twilight flats. Both images show the counterimage 
component to be extended, most probably due to blending with the lensing galaxy. The images are 
shown in Figure 4, where we have also attempted to subtract our best estimate of the contribution 
of the quasar images from the I-band image (see Section 3.2) to reveal the lensing galaxy. Table 1 
lists the broad-band magnitudes, flux densities and image-to-counterimage ratios. 
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2.3. Radio imaging 

We observed J 1004+ 1229 at 1.67 GHz with the VLBA'^ on 8-9 October 2001, using ah 10 
antennas. We used the compact radio source J1002+1216 as the phase reference (26' away, peak 
flux density 0.13 Jy), with a switching time of 5 minutes. The total dwell time on J1004+1229 was 
200 minutes. 

The data were recorded in eight 8 MHz bands in dual-polarization mode, for a total observing 
bandwidth of 32 MHz per polarization. Each 8 MHz band was subdivided into 16 x 0.5 MHz 

channels. The sampling time was 1 s. For the correlation phase centers, we used the JVAS position 
(J2000) of the phase reference source (10t'02™52!8457, +12?16™14'.'588; Browne et al. 1998), and 
the FIRST position of the target source. 

Calibration was performed with standard AIPS tasks. We used a fringe-fitting interval of 2 
minutes. Prior to imaging, we averaged the data in time into 6 s bins and in frequency into 1 MHz 
bins. This level of sampling was suflicient to reduce bandwidth and time-average smearing to less 
than 1% within 1''5 of the phase center. Imaging was also performed with AIPS. We created a 
2" X 2" image centered on the expected A-B midpoint. After component A was subtracted using 
the CLEAN algorithm, a significant peak (~ 9a) was apparent in the residual image. This peak is 
the brightest peak in the residual image, it does not belong to the sidelobe pattern of component A, 
and it is near the expected location of component B. We concluded that it represents component B 
of the lensed system. The final contour maps, shown in Figure 5, were produced after phase-only 
self-calibration with a 60 minute solution interval. 

When A and B are modeled as elliptical Gaussian components, their flux densities are 6.23 mJy 
and 0.99 mJy, respectively. The noise on the map is 0.11 mJy beam~^, and the the flux density 
ratio (A/B) is 6.3 ± 0.7. The overall flux density calibration is uncertain by 5%. Component A 
is 263.3 ± 1.0 mas west and 1517.2 it 1.7 mas south of component B. The position of A is (J2000) 
10l'04^'^24!8858, +12?29'.°22'.'313. The lowest surface-brightness contours of component A suggest 
the source is partially resolved. The elongation is in a direction nearly perpendicular to the A/B 
separation, which is suggestive of the tangential stretching that is characteristic of gravitational 
lensing. However, it is possible that the elongation is an artifact of residual gain errors rather 
than intrinsic source structure. Further evidence for extended emission comes from the difference 
between the 1.4 GHz flux density from FIRST of 12.3mJy, measured in an 5.6-arcscc beam, and 
the 7.2 mJy total flux measured at 1.67 GHz in the VLBA image. The source has a flux density at 
1.4 GHz of 11.8 mJy in the « 45-arcsec resolution NVSS survey (Condon et al. 1998), identical to 
the FIRST flux density to within the errors. Thus assuming the difference in flux density between 
the VLA and VLBA measurements is not due to variability, this suggests the VLBA has resolved 
out extended emission from the radio source, and also that the source is not signiflcantly larger 



''The Very Long Baseline Array (VLBA) is operated by the National Radio Astronomy Observatory, a facility of 
the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. 
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than the FIRST beamsize. 

3. Results and Interpretation 

3.1. The quasar spectrum 

Our ESI spectra of J1004+1229 show a comphcated absorption hne spectrum with an apparent 
edge at 10200A, which we ascribe to Mgll absorption at ^; = 2.65 (Figures 2 and 3). This interpre- 
tation is confirmed by the K-hand spectrum, which shows a strong, broad (3500 kms~^) emission 
hne at 2.40//m, corresponding to Hq: at z = 2.66 it 0.01. Most of the absorption features can be 
tentatively identified with low- ionization species; for example, Mgii 2796/2803 and Aim 1859 are 
clearly present. Much of the absorption between these wavelengths is due to mctastablc excited 
states of Fe II and perhaps Fe III (Hazard et al. 1987; Becker et al. 1997), hence its classification 
as an FeLoBAL. The absorbing systems which give rise to the spectral features seem to be quite 
complex, making the resemblance between J1004-I-1229 and FIRST 1556+3517 particularly striking 
(Figure 3). Nevertheless, there arc detailed differences in the spectra, for example the Fe I 2167 
line, apparently present in FIRST 1556+3517 is absent from J1004+1229. The blucshifts extend to 
« 10000 kms~^ even for the low ionization lines. J1004+1229A has very black absorption troughs, 
indicating a high covering factor for the absorbing clouds. 

3.2. Evidence for lensing 

The infrared and optical spectra, which indicate that both components A and B arc members of 
the rare FeLoBAL class of quasar, taken together with the similarity of the image-to-counterimage 
ratio in the near infrared and the radio indicate that components A and B are almost certainly 
lensed images of a single quasar. The optical spectra of A and B do, however, differ in detail, and 
in this section we discuss whether these differences are consistent with lensing. 

In a simple, isolated three-image lens system, the brightest component is farthest from the 
lensing galaxy, and the faintest (demagnified) component is closest. As third images are only rarely 
detected we modelled the system as containing two images, and assumed that the likely effects of 
lensing on the spectrum were blending of the weaker component with light from the lensing galaxy, 
and possibly differential reddening of the two images by dust in the lens. This picture is borne 
out by the photometry (Table 1): the image-to-counterimage ratio increases from the radio to K' 
and from K to J, as one would expect if the counterimage is more strongly reddened, but then 
decreases again in / as light from the lensing galaxy is blended into that from the counterimage. 

We thus attempted to decompose the spectrum of component B into a reddened version of 
component A plus the lensing galaxy. Fortunately the large number of strong features in the 
spectrum of the quasar makes this relatively easy to do. The unknown nature and redshift of the 



- 7- 



Ha 




2.0 2.2 2.4 

Wavelength / jxm 

Fig. 1.— K-hand spectrum of J1004+1229, obtained with NSFCAM at the IRTF, showing the Ha 
emission line in the two components A (lower) and B (upper). 



Band Mag/Flux A/B ratio 

V 22.9 ±0.1 ^4 

/ 19.7 ±0.1 4.1 ±0.5 

J 16.55 ±0.16 7.7 ±0.2 

K' 14.54 ± 0.09 6.8 ± 0.2 
1.4 GHz 12.3 mJy 

1.67GHz - 6.3 ± 0.7 

Table 1: Photometry and image ratios for the J1004+1229 system. Near-IR photometry is from 
the 2-MASS catalogue, optical photometry from our NOT observations and the radio flux density 
from FIRST. Image-to-counterimage ratios were determined from PSF fitting using DAOPHOT on 
our NOT and IRTF data in /, J and K . The signal-to-noise ratio in the y-band image was too 
low to permit a reliable measurement of the A/B flux ratio. 



Quasar redshift 2.65 ± 0.01 

Quasar Mb -28.5 + 2.51g(/Lt/3) - 2.51g(/i3 /3) 

1.4GHz radio luminosity 2 x 10^^ x (3//i)WHz-^sr-^ 

Lens redshift 0.95 ± 0.01 

Lens Mb ~ -21.1 

Lens (Tt, ^ 240kms~-^ 

Lens M/Lb « HHMq/Lq 

Table 2: Derived quantities for the J1004±1229 system. The lensing magnification of the quasar 
is denoted fi, and dust attenuation factor in rest- frame i?-band fs- The radio flux is assumed to 
originate entirely from the quasar. The mass-to-light ratio is quoted for Hq = 100/ikms~^Mpc~^ 
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Fig. 2. — The spectrum of components A and B of J1004+1229. The spectral resolution in these 
smoothed spectra is sA. The dotted lines, labelled in the upper panel, mark the redshifted wave- 
lengths of some of the stronger broad absorption features in the spectrum. 
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Fig. 3.— The spectrum of J1004+1229A compared to that of the FIRST FeLoBAL J1556+3517 
and the original FeLoBAL, Q0059-2735 (Hazard et al. 1987). The similarity of J1004+1229 to 
J1556+3517 is particularly strong. Prominent absorption features typically found in FeLoBALs 
are marked; the asterisk denotes a metastable state. 
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Fig. 4. — Images of J1004+1229. Top row, left to right: F, /, and the /-band image with our best 
estimate of the the two lensed quasar images subtracted to reveal the lensing galaxy. Bottom row, 
left to right: J and K' . Crosses on the /-band image mark the positions of the radio components 
detected with the VLB A (see Figure 5). 
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Fig. 5. — The 1.7 GHz VLBA image of J1004+1229. The bottom left panel shows the whole map, 
with the bottom right and top left panels showing close-ups of components A and B respectively. 
The dirty beam is shown in the top-right panel. 
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lensing galaxy mean that a large parameter space would need to be covered in a full analysis; we 
have therefore decided to restrict ourselves to a fairly course grid of reddenings to simplify the 
problem. 

We explored reddening values of E(B — V) from to 0.5 in steps of 0.1, assuming a ratio of 
selective to total extinction, Ry = 3.1 and a Milky Way extinction law (CardcUi, Clayton &: Mathis 
1989) (Figure 6). Initially we assumed a lens redshift zl = 1 at which we applied the extinction. 
We fixed the true image :counterimage ratio by estimating the excess reddening of component B at 
iC-band and then divided the K-hand image:counterimage ratio by this excess reddening factor. 

An E[B — V) ^ 0.2 was found to best remove the quasar features from the spectrum of B, 
leaving a residual which is smooth apart from a break at 7800 A. We interpret this as a 4000 A break 
at 2; = 0.95. In the lower panel of Fig. 6 we show the residuals from subtracting normalized, 
reddened spectra of component A from that of component B for reddenings of E{B — V) = 0.1, 0.2 
and 0.3 at the lens redshift. The intrinsic flux ratio we obtain for the best-fit, namely 5.5 is a little 
lower than the 6.3 obtained from the VLBA data (though only by 1.2-cr). Inspection of Figure 6 
seems to indicate that the intrinsic flux ratio lies between 5 and 6.2 (corresponding to reddening 
by the lens between E{B — V) = 0.1 and 0.3), but spatially-resolved spectroscopy of the lens is 
required to make a more accurate estimate. 

As the resolution of our images is insufficient to show the lens resolved from component B, 
we adopted a simple singular isothermal sphere model with which to extract the rough physical 
parameters of the lens system. In this model, the image to counterimage ratio is given by: 



where 9s is the angle between the source position and the lens center, and a is the deflection angle 
induced by the lens (which for this model is independent of impact parameter). The bright and 
faint images are offset by a + and a — 6s from the lens center, respectively. Thus given the 
image:counterimage ratio of 5.5 and separation (2a) of 1.54 arcsec, we expect the fainter image to 
be offset by 0.24 arcsec from the lens, consistent with us failing to resolve them in our images. The 
total magnification in this model is 2a /9s = 2.9. 

We derive a velocity dispersion for the lens of ~ 240km s^^, a mass within the Einstein ra- 
dius of 3.0 X 1O^^M0, and a blue mass-to-light ratio of ^ llh in solar units (where Hq = 
100^kms~^Mpc~^), all of which are consistent with the lens being a massive elliptical galaxy. The 
optical magnitude estimated from the spectra and the I-band image is also consistent with this, 
and suggests that the galaxy has a stellar luminosity about 0.5 magnitudes brighter than at 
zero redshift. 
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4. Discussion 

4.1. The dust in the lensing galaxy 

The low but non-zero value for the excess reddening produced by the lensing galaxy is at the 
upper end of the range observed for local ellipticals (e.g. Goudfrooij &: de Jong 1995), and also those 
inferred for other quasars lensed by 0.3 z 1 elliptical galaxies (Falco ct al. 1999). Elliptical 
galaxies at z ~ 1, however, are more likely to have low- level star formation episodes occurring in 
their cores (Menanteau, Abraham & Ellis 2001), so it is possible that such galaxies would typically 
have more dust in their inner few kpc than is usually seen either in local ellipticals or in z ~ 0.5 
lensing galaxies. 

4.2. The intrinsic nature of the quasar 

LoBALs are in general redder than normal quasars, almost certainly because of reddening in 
the quasar environment. Estimating the reddening of this quasar is difficult because of the strong 
absorption features in the rest-frame UV. If, however, J1004+1229 is like lower redshift FeLoBALs 
there should be little line absorption longward of Mgii 2798, and we can make a very rough estimate 
of the reddening by comparing the observed optical to near-infrared spectral energy distribution 
(SED) with that expected from a normal FBQS quasar at this redshift. We assume that the effect 
of reddening in the lens is negligible compared to that in the quasar environment: we expect dust 
only in the inner few kpc of the elliptical, and, as discussed above, this has only a small effect on the 
fainter image; the brighter image, with an impact parameter ~ 10 kpc is unlikely to be significantly 
reddened by dust in the lens. In Figure 7 we plot a comparison of the SED of J1004-I-1229 with that 
of FBQS composites (Brotherton et al. 2001) reddened with the Small Magellanic Cloud (SMC) 
extinction law of Pei (1992), from which we estimate an Ay ~ 1 mag. This is similar to the 
Ay ~ 1.5 estimated for the FeLoBAL J1556-F3517 by Najita, Dey & Brotherton (2001) (see also 
Fig. 1). 

With a radio luminosity at 5GHz, LsgHz lO^^WHz^^sr^^ after correction for lensing, the 
quasar is in the radio-intermediate class when this is defined through radio luminosity. In terms of 
the radio-to-UV continuum ratio R* (Sramek &; Weedman 1980), R* ~ 280 before correction of the 
UV flux for reddening and BAL absorption. After correction by an assumed E{B — V) = 0.35 for 
reddening, and an estimated BAL absorption of a further factor of two, this results in an R* 10, 
again placing the quasar in the radio-intermediate class, as is typical for BALs in FIRST (Becker 
et al. 2000; 2001). 

At restframe optical wavelengths the extinction due to dust approximately matches the mag- 
nification due to gravitational lensing, and using this assumption we can use the width of the Ha 
line to estimate a black hole mass of a few times 10^ Mq, following the prescription used by Lacy et 
al. (2001). Even with this high black hole mass, the high intrinsic luminosity of the quasar requires 
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Fig. 6. — Decomposition of the Keck ESI spectrum. Top panel: raw spectrum of component B 
(middle spectrum) compared to the spectrum of A divided by the iC-band flux ratio of 6.8 (upper 
spectrum) and the spectrum of A reddened by E{B — V) = 0.2 at z = 0.95, divided by 5.5 (lower 
spectrum). Lower panel: residuals obtained by subtracting from the spectrum of component B 
spectra of component A reddened in the rest-frame of the lens and normalized to the K-hand flux 
ratio. From top to bottom: A reddened by E{B — V) = 0.1 and divided by 6.2; A reddened 
by E{B -V) = 0.2 and divided by 5.5, and A reddened by E{B - V) = 0.3 and divided by 5. 
Superposed in red on all three spectra is a 6Gyr old, IGyr burst galaxy model from the PEGASE 
library of Fioc & Rocca-Volmerange (1997) redshifted to z = 0.95, our best estimate of the lens 
redshift. 
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that it is accreting near the Eddington limit. The high black hole mass is consistent with its high 
radio luminosity for a BAL quasar. 



5. Lensing of red quasars 

The steep number counts of high luminosity quasars produces a large magnification bias, i.e. 
a tendency for lensed quasars to be over-represented in a flux-limited quasar survey (e.g. Turner, 
Ostriker & Gott 1984). This typically results in such samples containing ~ 1% of lensed quasars, 
compared to ~ 0.1% expected in the absence of magnification bias. As our 2MASS quasar sample is 
very bright, the magnification bias is even higher than usual, and, using the lensing model discussed 
below, we estimate that the expected lensing rate for z > 0.5 quasars in our sample is « 3%. Even 
so, the discovery that 2/13 of our z > 0.5 luminous red quasars are lensed is surprising. Although 
the numbers of quasars and lenses in our sample is small, the Poisson probability of finding two or 
more lensed quasars in our sample by chance is only 0.06. Why the lensed fraction is so high is not 
clear, but it is probably linked to the observation of Chartas (2000) who showed that the fraction 
of lensed quasars which are BALs is higher than that in normal quasar samples. Chartas argues 
that the fraction of attenuated BAL quasars in samples of gravitationally-lcnscd quasars can be 
explained as a result of the joint effects of lensing and attenuation on the selection of these samples. 
His analysis is strongly dependent on the relative amounts of attenuation by dust and magnification 
by lensing, in addition to the slope of the luminosity function. To see this, we follow Goodrich 
(1997) and Chartas (2000) by considering a quasar population with a luminosity function <p{L, z) 
studied in a narrow range of luminosity (Li, L2) at redshift z. We draw samples of lensed quasars 
magnified by a mean magnificaton factor M, and reddened quasars (including BALs), attenuated 
by a mean attenuation factor A. We also now define the probability of the quasar being lensed by 
pi^ and the probability of it being significantly reddened by Pr- The fraction of lensed quasars in 
the sample is 

„ ^Li/M,L2/M,z) , 
fi = T-TT — p T Pi = bipi 

where ^ is the integral luminosity function and bi is the usual lensing magnfication bias factor (e.g. 
Turner, Ostriker & Gott 1984). The fraction of red quasars in the sample is fr'- 

_ ^LiA,L2A,z) _ 

where we define hr as the "attenuation bias factor" , by analogy with hi . If the luminosity function 
is steep (i.e., equivalent to a power-law of index < — 1 over the luminosity ranges of interest), then 
6; > 1 and 6^ < 1- The fraction of lensed quasars which are reddened, fir is then: 

_ ^{LiA/M,L2A/M,z) 
^Li/M,L2/M,z) 

Hence if A M and the luminosity function is a single power-law, then fir = Pr/h ~ hrPr = fr-, i-e. 
the same as in the rest of the sample. It is, however, suggestive in this context that both the lensed 
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Fig. 7. — The spectrum of J1004+1229A (solid blue line) compared to the FBQS composite spec- 
trum (Brotherton et al. 2001) reddened by E{B - V) = 0.3 and E{B - V) = 0.4 (upper and 
lower red dotted lines respectively). Points from the 2-MASS J,H, K and our NOT V and /-band 
photometry of the system are shown as green stars. The optical and infrared spectra have been 
scaled to match the broad-band magnitudes in / and K respectively. 
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objects are amongst the brightest in K-hand of all our FIRST-2MASS quasar candidates, where 
bi should be highest. A steep cutoff in the luminosity function at high luminosities could have the 
effect of dramatically decreasing br (as there are very few high luminosity quasars to attenuate) 
whilst keeping bi moderate. No evidence for such a cutoff has been found, however (e.g. Wisotzki 
2000). 

We now make the discussion more quantitative, by considering hypothetical quasar surveys 
selected in the K-hand, with magnitude limits oi K = 15, typical of the deeper portions of 2MASS, 
and K = 18, which would be typical of future large-area near-infared surveys with 4-m class 
telescopes, and S-band selected samples at B = 18 and B = 21. We adopt the luminosity function 
of Boyle et al. (2000) ^, and assume an optical spectrum whose flux density fi, oc We take the 

lensing probability pi{z) from Kochanek (1993), and integrate over the magnification distributions 
of King &; Browne (1996). The distribution of attenuations in the quasar population is uncertain. 
The only sample for which this has been well-studied is the extremely radio-bright 3C sample, 
which, of course, may well not be typical of the quasar population as a whole, and contains no 
BALs. Simpson, Rawlings & Lacy (1999) find that, in a total sample of 32 z ~ 1 3C radio sources, 
28^^|% of quasars are significantly reddened with Ay values ranging from ~ 2 to ~ 15. Given the 
difficulty of defining a plausible attenuation distribution we have decided to investigate how the 
fraction of reddened quasars in both lensed and unlensed samples varies as a function of attenuation, 
using the SMC extinction law of Pei (1992). In Fig. 8, we plot the fraction of reddened quasars 
expected in magnitude limited samples as a function of attenuation in the rest-frame at z = 2 
in both lensed and unlensed samples, selected in i^'-band and S-band. The plots are normalized 
to unity for zero attenuation. Even for samples selected in the near-infrared, red quasars with 
rest-frame Ay ~ 1 are strongly selected against. This is particularly true in bright samples, where 
the slope of the luminosity function is steepest. Also, although the fraction of red quasars in lensed 
samples is higher than that in unlensed samples, it remains a factor of several lower than the true 
fraction. These numbers depend only on the magnitude limit of the quasar sample and the band 
to which that limit is applied, and are independent of other selection criteria which might further 
decrease the number of red quasars in a sample, e.g. a blue color selection criterion. Our model 
implies that the high fraction of red quasars seen in lensed samples is not likely to be entirely due 
to magnification bias effects, if the assumption that reddened and unreddened quasars are drawn 
at random from the same underlying population is correct. In the S-band, where quasar samples 
are typically selected, the selection effects are even more pronounced, and in a B 18 magnitude- 
limited sample (e.g. LBQS) very few z > 2 quasars reddened by Ay > 1 in the rest-frame would 
ever be found. 



*At x ~ 2 this luminosity function is consistent with that of Wisotzki (2000), which, although it is only applicable 
to the bright end, covers the range of magnitudes which includes our very luminous quasars. 
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Fig. 8. — The fraction of red quasars in magnitude-limited samples at z = 2 as a function of rest- 
frame extinction. Left: the solid and dotted curves correspond to normal and lensed quasar samples 
in = 15 (red) and B = 18 (blue) samples. Right: the same for a sample with a. K = 18 (red) 
and a B = 21 (blue) magnitude limit. For comparison, J1004-I-1229 has Ay ~ 1 and J0134-0931 
(lensed) and J0738-h2927 (not lensed; Gregg et al. 2001) have Ay ~ 2. 

An alternative explanation for the high fraction of red lensed quasars is that red quasars are 
a distinct population, as might be the case if redness is a phase in a quasar's life cycle rather than 
a purely geometrical effect. This is consistent with the interpretation of BAL QSOs, especially the 
FeLoBALs, as being an early or heavily obscured phase in the lives of quasars (Voit et al. 1993; 
Egami et al. 1996; Becker et al. 1997, 2000; Gregg et al. 2000). If the luminosity function for 
red quasars is significantly steeper than for normal quasars, then this could lead to a higher than 
normal magnification bias in samples of red quasars, and explain our seemingly large lensing rate, 
as well as the tendency for the fraction of BAL QSOs in lensed samples to be higher. A physical 
basis for this might be that the luminosity function of quasars accreting at close to the Eddington 
rate will resemble the black hole mass function, which has a steep slope at high masses caused by 
the cutoff in galaxy masses produced by the Schechter function (e.g. Salucci et al. 1999). The only 
way for a geometrical model to be consistent with this result is if the range of angles over which 
the quasar is seen as significantly reddened increases as quasar luminosity decreases. 



6. Conclusions 

The comparison of FIRST and 2MASS has allowed us to discover quasars reddened by dust, 
both intrinsic and intervening. J1004-I-1229 is almost entirely reddened in the local environment 
of the quasar by a combination of dust and strong low ionization absorption features; the lens 
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contributing only a small amount of reddening, and probably only to the fainter lensed image. 

Even after taking into account magnification by lensing, the quasar is intrinsically very lumi- 
nous and contains a supermassive black hole, probably accreting at close to the Eddington rate. 
The apparently high frequency of lensing among our red quasars could be explained in a model in 
which red quasars have an intrinsically steeper luminosity function than the normal quasar pop- 
ulation, although if the geometry of the reddening material is such that the range of angles in 
which reddening is seen increases as quasar luminosity decreases, a geometrical explanation could 
also work. In either case, the extent to which high-redshift, dust-reddened quasars in both lensed 
and unlensed are selected against in bright quasar surveys is very great, and suggests that a large 
fraction of quasars are missed from surveys with optical or near-infrared selection. In particular it 
suggests that FeLoBALs, selected against both because of dust reddening and extreme rest-frame 
UV absorption, may be much more common than their rarity in magnitude-limited samples would 
suggest. When we complete our FIRST/2MASS survey, and analyse the results in the context 
of our modelling we should obtain a good estimate of the true fraction of red quasars in quasar 
samples of much more moderate radio luminosity than those in the 3C survey. 

The lensed nature of J1004-I-1229 gives us a unique opportunity to study the host galaxy of 
a high-redshift BAL through deep imaging with NICMOS on HST, or with ground-based adaptive 
optics. If red quasars and BALs are indeed a distinct and early phase in quasar evolution, study 
of the host galaxy of J1004-I-1229 may help reveal the phenomena responsible for creating such 
objects. 
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